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We study the dynamical breaking of chiral and color symmetries of dense quark matter in the 
ultrastatic hyperbolic spacetime R Cg> H 3 in the framework of an extended Nambu-Jona-Lasinio 
model. On the basis of analytical expressions for chiral and color condensates as functions of 
curvature and temperature, the phenomenon of dimensional reduction and gravitational catalysis 
0^ ' of symmetry breaking in strong gravitational field is demonstrated in the regime of weak coupling 

constants. In the case of strong couplings it is shown that curvature leads to small corrections 
to the flat-space values of condensate and thus enhances the symmetry breaking effects. Finally, 
using numerical calculations phase transitions under the influence of chemical potential and negative 
(3JT). curvature are considered and the phase portrait of the system is constructed. 
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I. INTRODUCTION 

Dynamical breaking of chiral and color symmetries has been successfully studied within field theories of the Nambu- 
Jona-Lasinio (NJL) type with four-fermion interactions They are also quite useful in describing the physics of 
Qh| light mesons (see e.g. @, 0, Q and references therein) and diquarks 

The possibility for the existence of the color superconductive(CSC) phase with a nonzero colored diquark condensate 
i ' was proposed both in the region of high baryon densities @, H Hi an d moderate densities [l(| HH, EH EE EH ■ In the 
framework of NJL models the CSC phase formation has generally been considered as a dynamical competition between 
diquark (qq) and usual quark-antiquark condensation (qq) . Special attention has been paid to the catalyzing influence 
of external fields on chiral symmetry breaking (xSB) in the regime of weak coupli ng | 15l. HH (constant magnetic 
field), [H, Eil (chromomagnetic fields) and on the condensation of diquarks [2(J l2l| (chromomagnetic fields). In 
particular, it was demonstrated that in a strong field the considered symmetry is dynamically broken for an arbitrary 
. weak attraction between quarks. The physical explanation for this is that the effect of dynamical symmetry breaking 
is accompanied by an effective lowering of dimensionality in strong fields, where the number of reduced units of 
dimensions depends on the concrete type of the background field. Dynamical symmetry breaking in a magnetic field 
' in spacetimes of dimension higher than four was also considered in [221 ] . 

For cosmological and astrophysical applications, it is also interesting to study the influence of spacetime curvature 
J* ■ on symmetry breaking. One of the ways to account for the effects of gravity is to use the adiabatic expansion of 
Green's functions in the vicinity of a fixed point in powers of small curvature (see, for example, the review [23|). 
However, since second-order phase transitions take place in the infrared region, the considered processes may become 
sensitive to the global structure of spacetime, and then one needs exact expressions for the propagators in the curved 
. i spacetime. 

In particular, exact solution can be found for spacetime with high symmetry. A variety of examples of x^B in 
different symmetric spaces both at zero and finite temperature has been considered in the literature (see e.g. f23|). 
One of the well-known examples of spaces with constant positive curvature is the Einstein universe of the form 
R <g> 5* 3 . xSB at finite temperature and chemical potential in the static Einstein universe was recently considered in 
|24l ]. Further investigations of quark matter in this gravitational background concerning, in particular, diquark and 
pion condensation were performed in [25, 26]. There, the positive curvature was shown to lead to the restoration of 
broken symmetries, thus acting in a similar way as the temperature. Another symmetric space, where exact solution 
may be found, is the ultrastatic hyperbolic spacetime R <g> H 3 . It was noted that in spaces with negative curvature 
the chiral symmetry is broken even at weak coupling constant (see e.g. [23]). The detailed analysis of the heat kernel 
in hyperbolic space showed that the physical reason for this phenomenon is the effective dimensional reduction for 
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fermions in the infrared region (2?j]. Recently, the combined influence of gravitational and magnetic fields on %SB 
in the special case of the 2D space of negative curvature, i.e., on the Lobachevsky plane, was considered in [28|]. As 
discussed by these authors, the study of the effects of surface curvature may be important for some condensed matter 
systems concerning, in particular, the quantum Hall effect in graphene [29l [30|. 

The aim of this paper is to study the effects of dynamical breaking of chiral and/or color symmetries in dense quark 
matter under the influence of negative curvature of ultrastatic hyperbolic spacetime R £g> H 3 . In the framework of 
an extended Nambu-Jona-Lasinio model, including (qq)- and (gq^-interactions, an exact in curvature expression for 
the thermodynamic potential is derived, which contains all the necessary information about the condensates. Basing 
upon the analytical solutions of gap equations for quark and diquark condensates we show that in strong gravitational 
field there arises a gravitational catalysis of dynamical symmetry breaking and chiral and color symmetries may 
be simultaneously broken even for weakly interacting quarks. This situation resembles the influence of magnetic or 
chromomagnetic fields on symmetry breaking in flat case. We also consider the role of finite temperature and find 
that for any fixed value of curvature there exists a critical temperature at which chiral and color symmetries become 
restored. Moreover, using numerical calculations, phase transitions under the influence of a chemical potential are 
investigated and the phase portrait of the system at zero temperature is constructed. Finally, it is shown that in 
the strong coupling regime negative curvature enhances the values of condensates as compared to the flat case. Our 
analysis demonstrates that negative curvature acts on chiral and color condensates in a way similar to that of a 
magnetic field in flat spacetime. 



II. THE EXTENDED NJL MODEL IN CURVED SPACETIME 

Let us briefly remind the basic definitions necessary for the description of fermions in curved spacetime. In 4- 
dimensional curved spacetime with signature (+, — , — , — ), the line element is written as 

The gamma- matrices 7/J , metric and the vierbein e^, as well as the definitions of the spinor covariant derivative 
V '„ and spin connection are given by the following relations [3ll |32|: 

{7/*(», 7*0*0} = 2g tiv {x), {7a, %} = 2rj &t , r?. g = diag(l, -1, -1, -1), 

W = *S> g^(x) = e^(x)e^(x), y„(x) = e%x) 1& . (1) 

V„ = d tl + r iM r„ = i w* 6 ^, <r aS = ~[7a,7s], 

— -^ e - aX z hp [C\p l i - C p a m - C^p], Ca pm = eld[ p e^a- (2) 

Here, the index d refers to the flat tangent space defined by the vierbein at the spacetime point x, and the j a (a = 
0, 1, 2, 3) are the usual Dirac gamma- matrices of Minkowski spacetime. Moreover 75, is defined, as usual (see, e.g., 
[33, [3J] ) , i.e. to be the same as in flat spacetime and thus independent of spacetime variables. 
The extended NJL model which includes the (qq)- and ((^-interactions of colored up- and down-quarks can be 
used to describe the formation of the color superconducting phase. For the color group SU C (3) its Lagrangian takes 
the form 



Gi 



> 2 



C = q [i 7 "V„ + M7°] Q + (qq) + (qi-ffq) + ^ H^q] [tqee\ 5 q c ] . (3) 



G 



2 L*« -,-b„,5„l r,v;wb„,5. 



Here, N c = 3 is the number of colors, G\ and Gi are coupling constants (their particular values will be chosen in 
what follows), [i is the quark chemical potential, q c = Cq l ,q c = q C are charge-conjugated bispinors (t stands for 
the transposition operation). The charge conjugation operation is defined with the help of the operator C = ij 2r y°, 
where the flat-space matrices 7^ and 7^ are used (see, e.g., [3l|). The quark field q = q la is a doublet of flavors and 
triplet of colors with indices % = 1,2; a = 1,2, 3. Moreover, f = (r 1 , t 2 , t 3 ) denote Pauli matrices in the flavor space; 
(e) lk = e lk , (e b ) Q/3 = e al3b are the totally antisymmetric tensors in the flavor and color spaces, respectively. 

Next, by applying the usual bosonization procedure, we obtain the linearized version of the Lagrangian with 
collective boson fields a, tt and A, 

t = q [i 7 "V M + M7°] q - Q + 9 - tJt(^ + ^) - 7^-A* b A b - A* b [iq'Ce^^q] - A b [iqee^Cq*] . (A) 
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The Lagrangians ([3]) and (|4|) are equivalent, as can be seen by using the Euler-Lagrange equations for the bosonic 
fields, from which it follows that 



Gi . , b 5 Gi _ _ Gi_. 5 ^ 

A = iq Gee 7 q, cr = — —qq, tt = — —qij rq. 

3 



3 



(5) 



The fields a and 7? are color singlets, and is a color anti-triplet and flavor singlet. 

In what follows, it is convenient to consider the effective action for boson fields, which is expressed through the 
integral over quark fields 



exp{iS' cff (CT,7r,A 6 ,A* fc )} = N J [dq][dq]exp^i J <£xj~gt 



where 



S eff (a,^,A 6 ,A* 6 ) 



d x\J—g 



2Gi 



3A b A *6 

~G~2 



(G) 



(7) 



with S q standing for the quark contribution to the effective action. 

In the mean field approximation, the fields cr, tt, A b , A* 6 can be replaced by their ground state averages: (cr), (tt), 
(A b ) and (A* b ), respectively. Let us choose the following ground state of our model: 

(A 1 ) = (A 2 ) = (t?) = 0. 

If (a) 7^ 0, the chiral symmetry is broken dynamically, and if (A 3 ) 7^ 0, the color symmetry is broken. Evidently, 
this choice breaks the color symmetry down to the residual group SU C (2). (In the following we denote (cr), (A 3 ) 7^ 
by letters cr, A.) 

The quark contribution has the following form (for more details see [2ll ]) 



S q (a,A) 



-i InDet 



iV - cr + M7° 



InDet 



4|A| 2 + (-iV - cr + M7°)(«V - cr + w ) 



(8) 



Here, the first determinant is over spinor, flavor and coordinate spaces, and the second one is over the two-dimensional 
color space as well, and V = J^V^. 

Note that the effective potential of the model, the global minimum point of which will determine the quantities a 
and A, is given by S c s = — V e g J d 4 Xy/—g, where 



^eff = 



3cr 2 3AA* 



2G X 



G-2 



Vei 



Sn 



v = d xy/—g. 



(9) 



III. THERMODYNAMIC POTENTIAL 

In this work we will consider the ultrastatic spacetime R <g> H 3 with constant negative curvature. The metric is 
given by 



ds 2 = dt 2 - a 2 (d9 2 + sinh 2 9dil 2 



(10) 



where a is the radius of the hyperboloid which is related to the scalar curvature by the relation R = — % , and dil 2 is 
the metric on the two dimensional unit sphere. 

Let us next introduce the one-particle Hamiltonian H — —iaS/ + (77 , where a — 7°7- Using this operator the 
quark contribution to the effective action can be written in the following form (for more details see [25I]): 



So 



In Dct 



H 2 - (ft, + 



2 In Dct 



4\A\ 2 + (H-(,) 2 -p 2 



}■ 



(11) 



where po = ido, and we have summed over colors (the Det-operator now does not include color space). 

In order to calculate the effective action one needs to solve the equation for eigenfunctions of the Hamiltonian H 
in hyperbolic space 



(-ioV + cr7°)* A = A* A . 



(12) 
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Decomposing ty\ into upper tp\,l an d lower i/jx,2 components and further separating variables by settin g = 
/i,2(^)Xi,2j the solutions can be found in terms of hypergeometric functions (for details see e.g. [lEHl] an d |34T3g|). 
Let us now consider the diagonal element: 



tr(x| In H 2 - {po + p) 2 \x) 



dP E E ln [« A W) 2 " (PO + M) 2 ] 4im(^4im(^ 
£— ± slm 



(13) 



where I, m, p = aV A 2 — a 2 and s are the quantum numbers which characterize the spinor \- Here £ is the sign of 
the energy and we will now take X(p) > 0. Notice that the last sum does not depend on the point x due to the 
homogeneity of the space H . Therefore it can be evaluated at the origin = 0, when only the I = term survives 
[3~i| (for generality, we quote below the formulae referring to the N-dimensional space) : 



(14) 



where 



Q(p) - 



2^-2 r(/ + jv/2)r(zp + 1/2) 
"7^ r(*p + i + n/2) ' 



(15) 



Notice that the definition of the measure pn(p) (the density of states) coincides with that used in [3£ 
Thus we can calculate the contribution from the first term in (jlll) : 



InDet 



H 2 - (Po + pf 



Tr ln 



H — {po + p) =N f / d N xdt tv{x,t\\n[H 2 - (po + nY]\x,t) 



(16) 



where v is the spacetime volume of R (g> i? , and Nf = 2 is the number of flavours. 
The second term gives: 



ln Dct 



4\A\ 2 + (H - p) 2 - p 2 = vN f 2^ I / dp^ N (p)ln[A\A\ 2 + (X(p)+^) 2 -p 2 o}- (17) 



2n 



Thus the quark contribution to the effective potential reads 

c(po 



cff 



>2^ /^E / ^^(p){ln[(A(p)+^) 2 -p2]+21n[4|A| 2 + (A(p)+^) 2 -^]}. (18) 



In the case of finite temperature T = 1//3 > 0, the following substitutions should be made: 



dpo , s _ _* \ 

9^ rt2^'"-'' 



2?r 



2tt 



1 



n = 0,±1 ± 2, 



where u„ is the Matsubara frequency. Then the quark contribution to the effective potential (|18j) becomes the 
thermodynamic potential fl q : 



+ OC 



20 



9.,, = --—^ E E / dpp N (p){\n[(\(p)+Zp) 2 +u J l]+2ln[4\A\ 2 + (\(p)+Zp) 2 +u J l]}. (19) 



=-oo £=± 

Summing over the Matsubara frequencies we obtain the thermodynamic potential 

v 2 |A|2 



} 



fi(<r,A) - iVc (^" + ^-) -2^ i iV / (iV c -2)^ 00 dp MJ v(p){A(p)+Tln(l + e -»)-^) + Tin (l + e -«*(/0+/0) 
-2^N f l°° dpfi N (p)h(X(p) - m) 2 + 4|A|2 + V(A(p) + m) 2 + 4|AP + (20) 

+ 2T111 (l + e -/3V(>(P)-^) 2 +4|A|^ + 2Tln ^ + e -/3V(A(p)+p) 2 +4|A|^ ^ 
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The spectrum of the Dirac operator A depends only on one dimensionless parameter p. Instead of p one can introduce 
the quantity with the dimension of momentum p = p/a. Then the spectrum of the Dirac operator may be written as 



A(p) = Vp 2 + ° 2 = E p . (21) 
In the case of the 3D space, N — 3, the density of states is 

M3(P) = £±i (22) 
Thus the thermodynamic potential in the i?(g> H 3 spacetime becomes 1 : 

4, J™ dp (p 2 + -^){E P + Tin (l + e H»(*»-/0) + Tin (l + e^^) 
+ J(E P pf + 4|A| 2 + J(E P + pf + 4| AP (23) 







2 /• 




Lr2 / 





2Tln ( 1 + e -^(^~M) 2 +4|Ap\ + 2T i n (i + e -l3V(E P +^+MA\A } 



It should be noted that the thermodynamic potential is divergent at large p. Therefore we must use some regular- 
ization procedure. Since we interpret p as the module of the momentum, the easiest way to regularize the divergent 
integral is to introduce the momentum cutoff, p < A. Then the regularized potential looks as follows: 

n^(a, A) = 3 (j^- + ^ - A £ dp (p 2 + + Thl ^ + e ~^ P -^ +Tln ^ + e -P(E P+ ») 



+ ^(E p -p)2+4\A\2 + ^(E p + pr+4\A\2 (24) 
+ 2Tln (l + e-V(^^) 2 +4|A|^ + 2Tln A + e -/3 v / (^+M) 2 +4|A|^ ^ 

where we used the expression for the scalar curvature R = — %. This formula consists of two parts arising from the 
two terms in the first round bracket of the integrand: the first part is the same as in (3+l)D Minkowski spacetime, 
while the second one, which is linear in curvature, corresponds to the contribution from (1+1)D spacetime. Hence, 
when the contribution of the second term dominates over the first one, dimensional reduction by two units takes place. 

The values of condensates a and |A| correspond to the point of global minimum of the regularized thermodynamic 
potential and are determined as solutions of the gap equations 

-IT = 0> ejAT = °" (25) 

In the following sections we will consider the behavior of the condensates as functions of curvature, temperature 
and chemical potential. 



IV. ANALYTICAL SOLUTIONS 



A. Chiral condensate 



Let us first consider the case when chiral symmetry is broken while the color symmetry remains unbroken (a ^ 
and A = 0). Since the quark condensate appears even in the vacuum, we put for simplicity p = and T = 0. 
Then one can obtain the following expression for the effective potential 

A 4 3t 2 3 \r\ n A 2 7? 

VcS (a) = -v (x), Vo ( x ) = — --F( x )-U-G(x), x = J: g = -jd, r = ± (26) 



Note that after introducing the quantity with dimension of momentum p = p/a the spectrum of the Dirac operator (12 1 1 formally 
coincides with the usual dispersion relation in flat Minkowski spacetime and does not depend on curvature. At the same time, the 
density of states 1221 1 in the hyperbolic space differs from the usual measure of integration over momentum in Minkowski spacetime by 
an additional term depending on the curvature (in fact, this is the only curvature dependent part of the thermodynamic potential). To 
obtain the correct measure of integration over the continuum spectrum one needs to use the properly normalized cigcnfunctions in 1141 . 
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where 



F(x) = (2 + x 2 )y/T^-x 4 ]n^-^-^-, G(x) = VT+~^ + x 2 In 1 + 

X X 

The gap equation for the condensate a reads 



^ = ^TTx^+( l 4-x 2 )ln- 
g 12 



(27) 



(28) 



An analytical solution of this equation can be obtained only for small x -C 1 (<r <C A). Expanding the right hand side 
of (12"51) in x we have 



|r| 
12 



x 2 ) In ■ 



(29) 



We will consider three different cases: a) subcritical g, g < g c = 1, where g c is the critical constant in flat four 
dimensional spacetime; b) near critical g, when g — > g c — and c) overcritical g, g > 1. 

In the case of subcritical g, a nontrivial solution of the gap equation (f2"9"|) exists only if x 2 < yjj-. In the strong 

curvature limit a 2 <C the chiral condensate is given by 



do = 2Aexp 



12 




= 2A exp 




e-o: 









12tt 2 (1 - g) 



\R\Gi 



(30) 



One can distinguish two sub-cases in which the last expression is consistent with the above made assumptions: 

i) g < 1, r <C 1 or ii) g <C l, r ~ 1. The case ii) shows that the gap equation has a nontrivial solution even at arbitrary 

weak coupling constant. 

The expression (f30|) looks very similar to the chiral condensate in the two-dimensional Gross-Neveu model. After 
excluding the coupling constant from the effective potential by using the gap equation (j29|) . we obtain 



V eS (a) 



\RW 2 

16tt 2 



lii - 



l 



(31) 



and this is indeed (up to a dimensional factor) the effective potential of the Gross-Neveu model. Hence, we con- 
clude that in the case of subcritical coupling the strong gravitational field of hyperbolic space leads to the effective 
dimensional reduction from (3+1) to (1+1) (see also [27j|). 

One should also note that the non-analytical dependence of the chiral condensate on curvature in the exponent of 
(|30f looks quite similar to that in a magnetic field 



Co = 



eB 



■ exp 



27T 2 (l- g ) 



\eB\Gi 



(32) 



but the pre-exponential factors differ (see for example the review 37]). This fact demonstrates that in the case of 
hyperbolic space the negative curvature plays an analogous catalyzing role as the magnetic field does in flat space, 
where it gives rise to the catalysis of x^B even at arbitrary weak attraction between quarks. 
In the near critical regime g — > g c — the chiral condensate is just the constant 




R\ 
12 ' 



(33) 



where the curvature must be small <C A 2 . 



■yj <C m 2 we obtain 



In the overcritical regime g > 1, a nontrivial solution of (|29[) exists only if x 2 > In the weak curvature limit 

\R\ . , /' if 



(To = m 1 



24m 2 



O 



(34) 



where m is the solution of the gap equation at R = 0. It is seen that in this case the curvature leads to small analytical 
corrections to the flat-space value of chiral condensate. 
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Next, let us consider the influence of finite temperature on the behavior of the chiral condensate. The temperature 
dependent contribution to the thermodynamic potential looks as follows 



1 2 f c 



,/,.|i> 2 + S), IM i . , i, 



Expanding the logarithm into a series and performing the integration over momentum, we obtain 

v/ 7T 2 \ 24 ^ n ^ n nBa I 

\ n—l n—l / 

where K v (x) is the Macdonald function (modified Besscl function). Then the gap equation at finite temperature reads 



1 fT, — 2 , J r l 2m 1 + \/l + a; 2 |r| 2 

- = VI + x 2 + — -x- In — h -x"{h -h), (35) 

g 12 x 12 

where 

oo oo 

h(l3a) = -2j2i-l) n Ko(nf3a), I 3 (f3a) = -2 ^(-l)"^^). 

n— 1 n—l 

Here we consider only the most interesting case of subcritical coupling and strong gravitational field. Thus, as 
previously at g < 1 and a 2 -yj, we obtain for the chiral condensate 



a (T) = 2Aexp 



I2n\l-g) 
\R\Gi 



h(Pa (T)) 



a (0)exp[-Ji(/J<7 (r))]. (36) 



The function /i(a;) is positive and monotonically decreases which means that temperature leads to the restoration 
of broken symmetry. The critical temperature T c is defined by the condition ao(T c ) = which gives the BCS-like 
relation 

T c = Tr- 1 e c a o (0) ^ 0, 57^(0), (37) 



where cto(0) is given by (|30|) . It should be mentioned that the temperature dependence of the chiral condensate (|36|) 
is the same as in a constant chromomagnetic field [2(| . 

This brief analysis shows that in the case of subcritical coupling the negative curvature leads to the catalysis of 
%SB while in the overcritical regime it leads to the enhancement of the chiral condensate. One might expect the 
same behaviour in the case of a color condensate. The case of overcritical coupling will also be considered in the next 
section using numerical methods. 

B. Mixed phase 

Let us now consider a more general situation, when both condensates can take nonzero values. What concerns 
color symmetry breaking in flat spacetime, it is well known that diquark pairing and color superconductivity arise 
(for large quark number densities and, correspondingly, at large chemical potential) due to an instability of the Fermi 
surface so that the color superconducting state is energetically more preferable. For studying the influence of gravity 
on the formation of a color condensate in its most "pure form", we find it convenient, contrary to the flat case, in the 
following to take /x = 0. Then the effective potential at zero temperature can be written in the form 



A 4 , 3A o o 1 



r 



V cS (a, A) = '—v ( Xl y), v (x) = ^-x 2 + By 2 - - (F(x) + 2F(z)) - £ (G(x) + 2G(z)) , (38) 
where 

a 2 I A I m * I 2 I 7 A 1 ^ o 37f2 R for^ 
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and the functions F and G are the same as in (|27l 
equations obtained from 



The nontrivial solutions for condensates satisfy the following gap 



where 



3A = H(x) + 2H(z) 



B = H(z) 



H(x) = V1 + . 



12 



x 1 In 



12 



(K(x)+2K(z)), 



K(z), 



i + VT+- 



K(x) = In 



1 + VT+: 



X X 

Substituting (|4T|) in (jiO)) . we obtain a separate equation for the chiral condensate 



(40) 
(41) 

(42) 



3A-2B = H(x) + j^K(x) 



(43) 



which formally coincides with (|28f in the previous Subsection A, when 1/g becomes replaced by (3 A — 2B). As in the 
previous case, we will solve equations (|43|) and (j4Tjl in the limit of small condensates x <C 1 and y <C 1 and consider 
only the most interesting case of strong curvature x 2 <C yjj-, y 2 <C yjj- and subcritical couplings A > 1, B > 1, where 
the effect of gravitational catalysis is more clearly seen. In this limit we simply have H (x) = 1 and K(x) = In (2/x). 
The solutions of the gap equations are as follows 



a 2 = 4A exp 



24 

-- n (3A-2B -1) 



™*o = cr o+4|Ao| 2 = 4A 2 exp 



24 
r 



(44) 
(45) 



where the coupling constants must satisfy the inequalities: 3 A — 2B > 1 and B > 1. From this, the color condensate 
follows 



4|A | 2 =ml \ 1-exp 



24 

— (3A-3B) 



(46) 



which exists only if A > B (i.e. G2 > §Gi). Hence, both condensates may exist simultaneously in the region 
A > B > 1 (in this region the inequality 3 A — 2B = A + 2(A — B) > 1 holds automatically). 
Excluding the coupling constants from the effective potential one obtains 



V e ff(a, A) = 



_\R]_ 

48tt 2 



In rr 



1 



2m In 



1 



(47) 



The other possible stationary point of the effective potential is a = and A ^ 0. This type of solution may be 
obtained from the previous one by letting ctq = in (|45|) so that the combined condensate m 2 Q reduces to the tilded 
expression 4|Ao| 2 . 

The phase structure of the model is defined by the global minimum of the effective potential. There are four types 
of stationary points of the effective potential: (0,0), (?o,0), (0, A ) and (cto, A ). Let us denote the corresponding 
values of the effective potential at these points by v\ , V2 , v% and V4 . The effective potential is normalized in such a 
way that Vi = 0. The other values are 



\R 

48 



v ± = — "ttt Wo +1ml 



«3 = -^|A | 2 , 
6 



'-'2 



16 ' 



(48) 



where ao is given by (|30p . First of all we see that the minimum v± = of the symmetry case is higher than the other 
minima. This means in contrast to the flat case that for subcritical couplings the symmetric phase in hyperbolic 
space is unstable under the formation of different condensates and symmetry breaking. Since m 2 = 4|Ao| 2 , it is 
easily seen that V4 < V3 in the region A > B > 1. The mixed phase is the true vacuum of our model if V4 < v 2 
or a 2 + 2m 2 > 3cr 2 . Dividing both sides of the last inequality by a 2 and using (|3"0|) , (|4"4")) and we obtain the 
following condition 



exp 



24 

^(B-A) 



2 exp 



24 

rM-B) 



> 3. 
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Introducing a new variable u = exp j^(B — A) , < u < 1 for A > B > 1, we see that the above condition is just the 

simple cubic inequality u 3 + 2 > 3u, which is automatically satisfied for < u < 1. Thus, the condition A > B > 1 is 
sufficient for the mixed phase to be the true vacuum state of our model. 

In the opposite case B > A > f the chiral and color condensates cannot exist simultaneously and we need to 
compare vi and vz- The inequality vj, < vi leads to A — B > jj^ In | which contradicts our previous assumption. 
Therefore in the region B > A > f only the phase with broken chiral symmetry occurs, and the chiral condensate is 
described by ([50)1 . 

Let us also briefly consider the influence of finite temperature on the mixed phase. Proceeding in the same manner 
as in Subsection A, we obtain the following temperature dependence of condensates 

a 2 {T) = a 2 (0)exph2/ 1 ( / 3 ( ro(T))], (49) 
ml (T) = m2o(0)oq)[-2/iC9m.o(T))] ) (50) 

where c 2 (0) and m 2 (0) are given by (|44|) and (|45|) . respectively. In analogy with (l37|) . the critical temperatures for 
chiral and color condensates are 

T c ff =7T- 1 e c cro(0), T C A = ^ 1 e c m,o(0). (51) 

Since the inclusion of temperature results in a smooth diminution of condensates, the relation between couplings which 
determines the phase structure of the model should not change. Again both condensates may exist simultaneously if 
A > B which implies the relation T A > T£, Therefore at T < T° the mixed phase is realized with chiral and color 
condensate of the form 

|A (T)| 2 = i{myT)-a 2 (T)}. (52) 

Obviously, at Tjf < T < T^, chiral symmetry is restored, while the color symmetry remains broken 

|A (T)| 2 = |Ao(T)| 2 = ^^. (53) 

Finally, at T > T C A both symmetries become restored. 

The phase portraits of the system are presented in FigQ]for fixed values of coupling constants A and B, for A < B 
(left picture) and for A > B (right picture). The critical curves which separate different phases are described by 
corresponding formulas (|3"T|) , (|5T) for critical temperatures. 




FIG. 1: The phase portraits of the system showing critical curves as functions of \r\ for fixed values of couplings A and B: for 
A < B (left) and for A > B (right). 



In fact the solutions (|4"9"| , (f5"01) of the gap equations give only an implicit dependence of the condensates Co (T) and 
m*o(T) on temperature. Due to the complexity of the function I\{x), no analytical solutions for condensates can 
be found. However, equations (|4*9"1) , (fSTJ)) can be easily solved by using numerical methods, for example, an iterative 
procedure. When the temperature is fixed, the result of such procedure will depend only on the first step of iteration. 
The most obvious choice is to take the condensates at zero temperature as a first approximation. The schematic 
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FIG. 2: Condensates ao, Ao and m*o as functions of temperature inside the mixed phase. (Henceforth we will omit the 
subscript of condensates on the plots.) 

picture of condensates inside the mixed phase is shown in Figj2]for fixed values of <7o(0) and Ao(0). It is interesting to 
note that the condensates (To and m^o have the usual form of the BCS theory of superconductivity and decrease with 
temperature, while the color condensate Ao slightly grows at T < T° and has a maximum around T° . We have not 
specified the values of condensates on these curves since we have no phenomenologically motivated choice of coupling 
constants and curvature. Fig|2] demonstrates only qualitatively the behaviour of condensates as functions of the 
temperature at fixed curvature \r\ and couplings A and B. Moreover, since the condensates at nonzero temperature 
depend on the curvature only through the condensates at zero temperature, the increasing curvature will produce the 
appropriate stretching of curves in Fig[2] 

V. PHASE TRANSITIONS 

In the previous section the phase structure of the model was extensively studied in the regime of subcritical couplings 
by formally considering arbitrary ratios of coupling constants. In particular, in order to investigate the influence of 
curvature on condensates in "pure form", the chemical potential was taken to be zero. Now, in this section, we turn 
to the more general case of phase transitions at finite chemical potential (quark number density). 

In the following we want to compare our results with the well established case of flat spacetime. Since the role of 
a growing temperature in the restoration of broken symmetries was already demonstrated, our considerations will be 
restricted here to the case of zero temperature. The investigation of the phase transitions at finite chemical potential 
beyond the limit of small condensates is difficult to perform analytically, and hence numerical methods will now be 
used. 

In order to be able to make comparison with the flat case, we consider the following "more physical" relation of 
couplings taken from the instanton-motivated NJL-model [ll[ (see also (2lj): 

G 2 = ?Gi- (54) 

o 

Let us now fix g — G\A? /it 2 = 1.4 in such a way that the chiral condensate at zero curvature, R = 0, in the vacuum 
is equal to the usual value of the constituent quark mass in flat space (ao = 350 MeV) . Here the cutoff parameter is 
taken to be A = 600 MeV. 

In terms of dimensionless couplings A and B (see (|5!)|)). relation (f54"| corresponds to B = 2 A, and one might think 
that it excludes the existence of a mixed phase in accordance with our previous discussion. However, we should 
stress that the inequality B > A, which guaranties the absence of a mixed phase, was obtained only for subcritical 
couplings. Here we have fixed g — 1.4 and thus A = 1/g < 1, and the overcritical regime is realized. The difference 
between these cases is that in the subcritical regime the main contribution to the values of condensates is given by 
the two-dimensional part of the effective potential, while in the overcritical regime, the condensates receive their 
contribution mainly from the flat four-dimensional part. Therefore our previous arguments are not applicable in this 
case. 



11 



It was observed earlier that in a flat four dimensional spacetime the relation (|54|) between coupling constants leads 
to the absence of a mixed phase in a wide range of parameters [Tlj . In the previous section, we have found that in 
the overcritical regime finite curvature gives only small corrections to the flat-space value of chiral condensate (see 
(|34[1 ). Using numerical calculations, we will now analogously show that in a wide range of the values of curvature its 
contribution to condensates is small in comparison with their values in flat case. It is clear that such corrections can 
not change the phase structure of the model obtained in flat case in [llj | , and thus can not lead to formation of new 
phases. Therefore, in what follows, we assume that there is no mixed phase, where both condensates simultaneously 
take nonzero values. 

As is well-known, for increasing chemical potential there arises diquark pairing, whereas the chiral condensate 
becomes suppressed. 

The corresponding behavior of the chiral and color condensates as functions of the chemical potential at r = and 
|r| = 1 is shown in FigJ3] 
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FIG. 3: Condensates a and A as functions of \i (all quantities are given in units of GeV) at r = (left) and at \r\ = 1 (right). 

As is seen from Figj3J the critical chemical potential, at which the phase transition takes place, in flat space, r — 0, is 
fi c w 330 MeV, while at |r| = 1 \i c ~ 366 MeV. From this we can conclude that the critical line /j. c (|r|) which separates 
the two phases is a growing function of the curvature. The numerical study of the number density, which is the first 
derivative of the thermodynamical potential with respect to the chemical potential, shows that it is discontinuous at 
fi c . Thus, an increasing chemical potential leads to a first order phase transition. 

It should also be mentioned that for the limiting case of zero curvature our results for the value of the critical 
chemical potential [i c » 330 MeV and the maximum value of the color condensate A w 70 MeV (see FigJ3]) are in 
agreement with the results obtained in [38| for the same values of <7o and A (note that our value of the color condensate 
A is by definition two times less than condensate A in (38j). 

We can also examine the condensates a and A as functions of the absolute value of curvature \r\. The behavior of 
chiral and color condensates inside both phases is presented in FigSJ 



cr 



0.4 
0.3 
0.2 
0.1 



A 

0.08 r 



0.06 
0.04 : 
0.02 



0.2 



0.4 



0.6 



0.8 



1.0 



0.2 



0.4 



0.6 



0.8 



1.0 



FIG. 4: Condensates (in units of GeV) a at fi = 200 MeV (left) and A at fi = 400 MeV (right) as functions of |r|. 
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It is clear that with growing curvature the values of condensates slowly increase, which results in the enhancement 
of the symmetry breaking effects. As we have already mentioned above, in a wide range of curvature values the 
increment of the condensates is small with respect to their values at r = 0. Moreover, it is seen from Fig[4] that 
condensates grow with curvature almost linearly which is in agreement with our previous analytical consideration (see 

mi 

Finally, the r — fj,- phase portrait of our system is shown in FigO As it was already mentioned, the critical chemical 
potential, at which the first order phase transition takes place, slightly grows with increasing curvature. 
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FIG. 5: Phase portrait for the coupling relation G2 = f Gi. 



VI. SUMMARY 



In the framework of an extended NJL model we have investigated chiral and/or color symmetry breaking in dense 
quark matter under the influence of negative curvature of hyperbolic space by employing the thermodynamic (effective) 
potential of the system as a function of chiral and diquark condensates, temperature, curvature and chemical potential. 

Two different regimes of dynamical symmetry breaking have been studied. First, the regime of subcritical couplings, 
where the values of coupling constants are lower than their critical values in flat space, was considered. As is well 
known, in this situation in four dimensions, no symmetry breaking takes place, and gap equations have only trivial 
solutions. Unlike the flat case, in hyperbolic space, symmetry may, however, be broken even for an arbitrary small 
coupling constant. In the case of subcritical couplings we obtained expressions for the chiral and color condensates 
that depend non-analytically on curvature. Secondly, the regime of overcritical coupling constants was investigated, 
when the symmetry may be broken even in flat space. In this case it was shown that curvature leads to small analytical 
corrections which increase the flat-space values of condensates and thus enhances the symmetry breaking effects. 

It is interesting to note that in the subcritical regime of coupling constants the strong grav itational field of hyperbolic 
space serves as a catalyzing factor similar to the role of the magnetic field [HI, [To. Il7j| or chromomagnetic fields 
0, E, H3, HU in the effects of dynamical symmetry breaking. As we have explicitly demonstrated, the gravitational 
catalysis takes place for chiral and color condensates. It is also worth mentioning that the effect of gravitational 
catalysis is accompanied by a lowering of dimensions of the system. In the regime of subcritical couplings the 
solutions of gap equations look quite similar to those for the 2D Gross-Neveu model. Therefore, we concluded that 
the strong curvature of hyperbolic space leads to an effective dimensional reduction by two units (see also [27j| ) . 

As we have already mentioned, in the subcritical regime the negative curvature essentially changes the phase 
structure of the NJL model found in [3{| making the symmetric phase unstable under the formation of condensates, 
while in the overcritical case one expects only minor modifications to the phase structure obtained in flat space. 
Therefore we have extensively studied the phase structure in the subcritical regime. It is interesting to note that the 
overall phase structure depends only on the ratio of the inverse (dimensionless) coupling constants A and B, but not 
on the curvature. For subcritical couplings the mixed phase is realized only if A/B > 1, while if A/B < 1 only the 
chiral symmetry may be broken. The same critical ratio A/B = 1 was found in the flat space in the framework of the 
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random matrix model [4(| • It was argued by these authors that this relation is a consequence of global symmetries of 
the model. 

Moreover, the influence of finite temperature on phase structure was investigated. The phase portraits of the 
system for different relations between couplings were constructed. In particular, we have shown that for any fixed 
value of curvature there exists a critical temperature at which the phase transition takes place and symmetry becomes 
restored. 

Finally, using numerical calculations, we have investigated the phase transitions between xSB and CSC phases 
under the influence of chemical potential and curvature in the regime of overcritical couplings. It was demonstrated 
that similar to the flat case there arises a diquark pairing for increasing chemical potential, while the chiral condensate 
becomes suppressed. The phase portrait of the system at zero temperature was also constructed, and it was shown 
that the critical line /i c (|r|), separating the two different phases, is a growing function of curvature. The chiral and 
diquark condensates, a and A, acquire only small corrections due to curvature increasing the flat-space values of 
condensates and this leads to an enhancement of the symmetry breaking effects. 

The results of this paper, although describing a model situation with symmetry breaking in a hyperbolic space, 
may hopefully find further development in more realistic situations with phase transitions in quark matter under the 
influence of strong gravitational fields. 
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